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Parity violating electron scattering allows model independent measurements of neutron densities 
that are free from most strong interaction uncertainties. In this paper we present statistical error 
estimates for a variety of experiments. The neutron radius R n can be measured in several nuclei, 
as long as the nuclear excited states are not too low in energy. We present error estimates for R n 
measurements in 40 Ca, 48 Ca, 112 Sn, 120 Sn, 124 Sn, and 208 Pb. In general, we find that the smaller 
the nucleus, the easier the measurement. This is because smaller nuclei can be measured at higher 
momentum transfers where the parity violating asymmetry A pv is larger. Also in general, the more 
neutron rich the isotope, the easier the measurement, because neutron rich isotopes have larger 
weak charges and larger A pv . Measuring R n in 48 Ca appears very promising because it has a higher 
figure of merit than 208 Pb. In addition, i?„( 48 Ca) may be more easily related to two nucleon and 
three nucleon interactions, including very interesting three neutron forces, than J? ri ( 208 Pb). After 
measuring R n , one can constrain the surface thickness of the neutron density a n with a second 
measurement at somewhat higher momentum transfers. We present statistical error estimates for 
measuring a n in 48 Ca, 120 Sn, and 208 Pb. Again, we find that a n is easier to measure in smaller 
nuclei. 

PACS numbers: 21.10.Gv, 25.30.Bf, 24.80.+y 



I. INTRODUCTION 

Nuclear charge densities have been accurately mea- 
sured with elastic electron scattering and have become 
our picture of the atomic nucleus, see for example ref. 
PQ. These measurements have had an enormous impact. 
Unfortunately, neutron densities are not directly probed 
in electron scattering because the neutron is uncharged 
[55] . Our knowledge of neutron densities comes primarily 
from hadron scattering experiments involving for exam- 
ple pions [2] , protons [3H5] , or antiprotons [H [7j . How- 
ever, the interpretation of hadron scattering experiments 
is model dependent because of uncertainties in the strong 
interactions. Often there are uncertainties in reaction 
mechanism and in probe-nucleon interaction. Particular 
probes may have additional uncertainties. For example 
antiprotons may only probe the large radius tail of the 
neutron density, because of strong absorption. Therefore 
neutron densities, deduced from hadron scattering, could 
have significant strong interaction uncertainties. 

Parity violating electron scattering provides a model 
independent probe of neutron densities that is free from 
most strong interaction uncertainties. This is because 
the weak charge of a neutron is much larger than that 
of a proton [£f. Therefore the Z° boson, that carries 
the weak force, couples primarily to neutrons. In Born 
approximation, the parity violating asymmetry A pv , the 
fractional difference in cross sections for positive and neg- 
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ative helicity electrons, is proportional to the weak form 
factor. This is very close to the Fourier transform of the 
neutron density. Therefore the neutron density can be 
extracted from an electro- weak measurement [8]. How- 
ever, the Born approximation is not valid for a heavy 
nucleus and coulomb distortion effects must be included. 
These were calculated in ref. [5] by numerically solving 
the Dirac equation for an electron scattering in both the 
coulomb potential and a weak interaction axial vector 
potential. Many details of a practical parity violating 
experiment to measure neutron densities have been dis- 
cussed in a long paper jlOj. 

The Lead Radius Experiment (PREX) measures the 
parity violating asymmetry A pv for 1.05 GeV electrons 
scattering from 208 Pb at five degrees [llj . This measure- 
ment should be sensitive to the neutron r.m.s radius of 
208 Pb to 1% (±0.05 fm) and initially ran at Jefferson Lab- 
oratory in the spring of 2010. PREX should demonstrate 
this new technique to measure neutron densities. 

The neutron radius of 208 Pb, R n , has important im- 
plications for astrophysics. There is a strong correlation 
between R n and the pressure of neutron matter at den- 
sities near 0.1 fm -3 (about 2/3 of nuclear density) |12) . 
The larger the pressure of neutron matter, the more neu- 
trons are pushed out against surface tension in 208 Pb, 
and the larger is R n . Therefore measuring R n provides 
very interesting information on the equation of state — 
pressure as a function of density — of neutron matter. 
The equation of state is very important in astrophysics 
to determine the structure of neutron stars. 

Recently Hebeler et al. [13] use chiral perturbation 
theory to calculate the equation of state (EOS) of neutron 
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matter. Their EOS has important contributions from 
very interesting three neutron forces. We have detailed 
information on two nucleon forces from nucleon-nucleon 
scattering. Furthermore, we have some information on 
isospin 1/2 three nucleon forces from mass 3 nuclei ( 3 He, 
3 H) and proton-deuteron scattering. However, our ex- 
perimental information on three neutron forces is lim- 
ited. Uncertainties in three neutron forces dominate the 
error bars of Hebeler et al.'s EOS. From their EOS, and 
the above correlation of R n with pressure, they predict 
R n — R p — 0.17 ± 0.03 fm. Here R p is the known proton 
radius of 208 Pb. Therefore, measuring R n provides an 
important experimental check on fundamental neutron 
matter calculations, and measuring neutron densities can 
provide important constraints on three neutron forces. 

The correlation between R n and the radius of a neu- 
tron star rjvs is also very interesting [14]. Both depend 
on the EOS of neutron rich matter. In general, a larger 
R n implies a stiffer EOS, with a larger pressure, that will 
also suggest ?-jvs is larger. Note that this correlation is 
between objects that differ in size by 18 orders of mag- 
nitude from R n « 5.5 fm to r^vs ~ 10 km. However, 
R n depends on the EOS at nuclear density and below 
— some average of the surface and interior densities in 
208 Pb. While, rNs also depends on the EOS at higher- 
densities, since the central density of a neutron star is a 
few or more times nuclear density. Therefore, measuring 
both R n and r/vs determines the density dependence of 
the EOS. For example if R n is relatively large then the 
EOS is stiff at low densities. If at the same time r^s 
is small then the high density EOS is soft with a low 
pressure. This softening of the EOS with density could 
strongly suggest a phase transition to a high density ex- 
otic phase such as quark matter, strange matter, or a 
color superconductor. 

Indeed, if one can simultaneously measure masses and 
radii of neutron stars with a range of masses one can de- 
termine the full neutron matter EOS, pressure as a func- 
tion of density [15] . For example, very low mass neutron 
stars of say 1/3 of M Q in mass have central densities 
near nuclear density. Therefore for these stars, the in- 
formation in rjvs is closely related to the information in 
Rn |16j . However these low mass neutron stars may not 
exist since they are difficult to form. 

Recently there has been great progress in deducing r/vs 
(for say one solar mass and above stars) from X-ray ob- 
servations. Ozel et al. find rjvs is very small, below 10 
km from observations of X-ray bursts [17] . While Steiner 
et al. reanalyze the burst data and include observations 
of neutron stars in globular star clusters [T8] • They con- 
clude that rjvs is near 12 km. From the EOS implied by 
these observations they predict that R n — R p = 0.15±0.02 
fm in 208 Pb. Note that this prediction for a nuclear size 
is based on astronomical observations. If Steiner et al. is 
correct, the high density EOS is relatively stiff and there 
is little room for a significant softening due to a transi- 
tion to a high density exotic phase of QCD matter. This 
potentially very important Steiner et al. result can be 



tested by measuring R n . Finally, Demorest et al. have 
recently discovered a 1.97 M Q neutron star [19]. This im- 
portant observation rules out soft high density equations 
of state. 

The EOS of neutron matter is closely related to the 
symmetry energy S. This describes how the energy of 
symmetric nuclear matter rises as one goes away from 
equal numbers of neutrons and protons. To a good ap- 
proximation the energy of neutron matter is equal to the 
energy of symmetric nuclear matter plus the symmetry 
energy. The pressure of neutron matter depends on the 
derivative of the energy with respect to density. There- 
fore there is a strong correlation between R n and the den- 
sity dependence of the symmetry energy dS/dp, with p 
the baryon density. The symmetry energy can be probed 
in heavy ion collisions [20] , For example, dS/dp has been 
extracted from isospin diffusion data [21] where projectile 
and target nuclei with different proton to neutron ratios 
are collided and one observes how the neutron to proton 
ratio equilibrates. Often isospin diffusion data is ana- 
lyzed with transport models, since heavy ion collisions 
are so complicated. Measuring R n allows one to extract 
dS/dp in a way that is independent of the complex HI 
dynamics. 

The symmetry energy S helps determine the composi- 
tion of a neutron star. Neutron stars ( assuming a nu- 
cleon phase) are thought to be about 90 % neutrons and 
10 % protons and electrons. The proton fraction is de- 
termined by S. The larger 5* is the larger the proton 
fraction. If the proton fraction is larger than about 12 % 
than a neutron n near the Fermi surface can beta decay 
to a proton p and an electron e and conserve both energy 
and momentum. 

n^p + e + v e , (1) 

The would be followed by p + e — > n + v e , cooling the star 
by the radiation of a vv pair. This is called the direct 
URCA process. If R n is measured to be relatively large 
then dS/dp is large, at normal density, so that S is likely 
large at high density. This would strongly suggest that 
massive neutron stars will cool rapidly by direct URCA 

m 

Finally both the solid crust of a neutron star and the 
skin of a heavy nucleus are made of neutron rich mat- 
ter at similar subnuclear densities. If R n is large then 
the energy of neutron matter rises rapidly with density 
and this quickly favors a transition to the uniform liquid 
phase. Therefore, there is a strong correlation between 
R n and the transition density p t in a neutron star from 
low density solid crust to liquid core [23 . If R n is large 
then the transition density from solid crust to uniform 
liquid core is predicted to be low. Many neutron star 
observables depend on properties of the solid crust. In- 
deed, the neutron radius R n impacts a very large range 
of nuclear physics and astrophysics. 

In this paper, we explore possible parity violating ex- 
periments, in addition to PREX, to obtain further in- 
formation on neutron densities. First, in Section [H] we 
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discuss some general considerations for neutron density 
measurements. Next, in Section III we review the PREX 



experiment and the determination of its statistical error. 
Then in Section [TV] we present statistical error estimates 
for a number of other neutron density measurements. We 
conclude in Section [V] that several measurements are fea- 
sible. We hope that this paper will lead to discussion of 
the physics impacts of these additional parity violating 
measurements of neutron densities. 



II. GENERAL CONSIDERATIONS FOR 
PARITY VIOLATING NEUTRON DENSITY 
EXPERIMENTS 

In this section we discuss general considerations for 
neutron density experiments. These considerations are 
necessary to control statistical and systematic errors. 
First the parity violating asymmetry A pv is small, of or- 
der a part per million. Therefore one needs to accumulate 
large statistics. This likely requires intense beams, high 
electron polarization, thick targets, and large acceptance 
detectors. 

The need for large acceptance detectors is made more 
challenging by the further requirement of good energy 
resolution to separate out inelastically scattered elec- 
trons. This also requires a nuclear target with a relatively 
high first excited state. As a result neutron density mea- 
surements on deformed heavy nuclei such as Uranium or 
on odd A heavy nuclei may not be feasible. 

A thick target, of order 10% of a radiation length, 
is likely necessary. This target will produce significant 
losses from bremsstrahlung, and there will be significant 
beam heating. However, a high powered solid target may 
be simpler than a high powered liquid hydrogen target. 
Furthermore, the presence of a target backing or windows 
may not be a significant problem. This extra material 
may only slightly dilute the measured asymmetry. 

The desired momentum transfer q is likely set by the 
size of the nucleus. The statistical error, at fixed q, tends 
to be smallest for higher beam energies and smaller scat- 
tering angles. Thus the statistical error is often mini- 
mized by using the smallest practical scattering angle. 
We define a figure of merit FOM x for using parity vi- 
olation to measure a property x of the neutron density, 
such as the r.m.s. radius R n or the surface thickness a ra , 
as follows, 

FOM x = ^^1a pv (9, E)\ x (0, Ef. (2) 

Here the differential cross section is da/dtt, the parity 
violating asymmetry is A pv , and e x is the sensitivity of 
A pv to changes in the quantity x, 
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is maximized for large beam energy E and small labo- 
ratory scattering angle 9. Note that A pv and e x depend 
primarily on q. However they also depend somewhat on 
E and 6 separately because of Coulomb distortions. As 
we calculate in the next sections, the statistical error in 
determining a quantity x is proportional to (FOM I ) -1 '' 2 . 

In addition to statistical errors one has to control er- 
rors from several other sources including helicity corre- 
lated beam properties, normalization, and radiative cor- 
rections. Systematic errors from helicity correlated beam 
properties involve false asymmetries from beam proper- 
ties that change with electron spin. These errors may be 
most important for experiments measuring small asym- 
metries. In general the asymmetry grows with increasing 
momentum transfer q. Therefore, errors associated with 
helicity correlated beam properties may be most impor- 
tant for measurements at small q. 

Normalization errors arise from imperfect knowledge 
of the electron beam polarization. The Lead Radius ex- 
periment may require beam polarimetry good to one % 
and a future, even more precise, experiment would likely 
require sub one % polarimetry. Finally, radiative cor- 
rections involve processes with additional real or virtual 
photons. They can be divided into coulomb distortion 
corrections where the nucleus remains in its ground state 
and dispersion corrections where the nucleus is in excited 
intermediate states. Coulomb distortion corrections have 
been accurately calculated in ref. [§] ■ Dispersion correc- 
tions have been estimated in a few cases, see for example 
|24) . Radiative corrections will become increasingly im- 
portant as the precision of a measurement is increased. 



III. THE LEAD RADIUS EXPERIMENT (PREX) 

In this section we review the Lead Radius Experiment 
(PREX) and calculate its statistical error. This will pro- 
vide context for the discussion of additional measure- 
ments in Section Hvl PREX aims to measure the neutron 
r.m.s. radius R n of 208 Pb with sensitivity to 1% (±0.05 
fm). We first collect some model neutron densities for 
208 Pb that are based on both nonrelativistic and rela- 
tivistic mean field interactions. These calculations yield 
a range of neutron radii as listed in Table |T] and shown in 
Fig. [T] We modify the well known NL3 relativistic mean 
field interaction |25j by adding a nonlinear coupling A„ 
between the rho b M and omega mean fields as dis- 
cussed in ref. [23] . The interaction Lagarangian density 
Ant is, 



£ int = A v g 2 p b„ ■ Wg 2 v V^. 



(4) 



If e x is large, one does not have to measure A pv as accu- 
rately in order to accurately determine x. Often FOM x 



The rho coupling constant g 2 p is also modified as listed in 
Table |H} This procedure was described in ref. [23] and 
generates neutron densities with a large range of neutron 
radii. The other parameters of the NL3 interaction are 
unchanged (23j [25] . 

Given these neutron densities, we calculate the parity 
violating asymmetry A pv by solving the Dirac equation 
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FIG. 1: (Color on line) Neutron density of 208 Pb for the 
Skyrme interactions SI (black dotted curve), SLY4 (short 
green dashes) and relativistic mean field interactions FSUg- 
old (longer red dashes), NL3 (solid blue), and NL3m05 (black 
dot dashed curve). See Table [I] and text. 



GeV and a scattering angle of five degrees, Fig. [3] plots 
A pv versus R n for the densities of Table [T] This figure 
shows that there is a simple relation between R n and 
A pv so that one can deduce R n from a measured A 



This will be discussed further below. 
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SI, R„=5.49 fm 
SLY4, 5.62 fm 
FSU, 5.68 fm 
NL3, 5.74 fm 
NL3m05, 5.85 
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TABLE I: Model root mean square proton R p and neutron 
R n radii for 208 Pb. For NL3p06 and NL3m05 see text and 
Table E] 



Interaction 


Rp (fm) Rn (fm) 


Skyrme I [26] 


5.38 


5.49 


Skyrme III [S7J 


5.52 


5.65 


Skyrme SLY4 [28] 


5.46 


5.62 


FSUgold [29] 


5.47 


5.68 


NL3 [H] 


5.46 


5.74 


NL3p06 


5.51 


5.60 


NL3m05 


5.50 


5.85 



as discussed in ref . [5] . Figure [2] shows A pv versus lab- 
oratory scattering angle for an electron beam energy 
of 1.05 GeV. Note that we use the experimental 208 Pb 
charge density [1] for these calculations instead of the 
model charge densities. Figure [2] shows that A pv is very 
sensitive to R n for scattering angles near five degrees, 
and A pv decreases with increasing R n . In contrast, for 
scattering angles between seven and eight degrees, A pv 
is mostly insensitive to R n and actually increases very 
slightly with increasing R n . For a beam energy of 1.05 

TABLE II: Isovector interactions for modified NL3 [25] pa- 
rameter sets, see text and ref. |23j . 

Interaction A„ g 2 p 
NL3 79.6 

NL3p06 0.06 300. 
NL3m05 -0.05 42.0 



FIG. 2: (Color on line) Parity violating asymmetry A pv ver- 
sus laboratory scattering angle for elastic electron scattering 
from 208 Pb at 1.05 GeV. The different curves are for neu- 
tron densities calculated with the indicated interactions from 
Tabled] 



G — O 850 MeV, 6 deg. 
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FIG. 3: (Color on line) Parity violating asymmetry A pv versus 
neutron r.m.s. radius R n for the model neutron densities of 
Table |T] The red squares are for a beam energy of 1.05 GeV 
and a laboratory scattering angle of five degrees while the 
black circles are for an energy of 850 MeV and a scattering 
angle of six degrees. 



We calculate the sensitivity e# n of A pv to the neutron 
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radius R n as follows, 



d In A™ 

e R n = 



R n dA. 



d In R n Ap V dR n 



(5) 



We calculate dA pv /dR n by calculating the change in A pv 
when a model neutron density p n (r) is streched by a fac- 
tor A 1.01, 



1 r 

Pn(r) -> T3i°n(-r)- 



(6) 



Figure [4] shows versus scattering angle for beam 
energies E of 1.05 and 1.8 GeV. Near E = 1.05 GeV and 
8 = 5 degrees en n peaks near 3. This shows that a 3% 
measurement of A pv is sensitive to the neutron radius to 
1% (±0.05 fm). 
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FIG. 4: (Color on line) Sensitivity of the parity violating 
asymmetry A pv for 208 Pb to changes in the neutron radius 
En n = eu n = din A pv /d In R n versus scattering angle 8 at 
beam energies of 1.05 GeV (solid line) and 1.8 GeV (dashed 
line). The SLY4 neutron density was used. 



We now estimate statistical errors assuming the exper- 
imental parameters in Table |III| appropriate for Hall A 
at Jefferson Laboratory [IT] where the two High Reso- 
lution Spectrometers have solid angle acceptances near 
A£l = 0.0037 Sr each. We somewhat arbitrarily assume 
a beam current I of lOOpA. This could be slightly opti- 
mistic. It is easy to scale calculated count rates for other 
beam currents. 

We assume a target thickness of 10% of a radiation 
length, see Table |IV[ and a modest energy resolution of 
AE = 4 MeV. This target thickness will introduce signif- 
icant radiation losses that depend somewhat on energy 
[30j . For simplicity we multiply our rate by a radia- 
tion loss factor £ = 0.34 to describe electrons that are 
not detected because of energy loss in the target and or 
bremsstrahlung during scattering. Note that our modest 



TABLE III: Experimental parameters including beam current 
I, beam polarization P, detector solid angle AO,, number of 
arms N, energy resolution AE, and radiation loss factor £, 
see text. 

Parameter Value 



I 
P 

An 

N 
AE 

c 



100^M 
0.8 

0.0037 Sr 
2 

4 MeV 
0.34 



TABLE IV: Target thickness t assuming 10% of a radiation 
length. The number of target atoms per square cm is ptar- 

Element t (cm) p tar (cm~ 2 ) 



208p b 


0.05 


1.6 x 10 21 


120 Sn 


0.12 


4.4 x 10 21 


48 Ca 


1.04 


2.4 x 10 22 



energy resolution of 4 MeV will include a small contami- 
nation from the 2.6 MeV 3~ first excited state in 208 Pb. 
However this inelastic cross section is small at forward 
angles and one can make a crude estimate of A pv for this 
3" state [TU]. 

The total number of electrons detected A^otai during a 
running time T is 



JVtotal = IT p U r ^CAftiV, 

dil 



(7) 



assuming a high efficiency detector. Here p tar is the num- 
ber density of the target in nuclei per unit area, see Table 
IV) This allows one to be sensitive to R n with a statistical 
error AR n , 



ARn 
Rn 



= ( iVtotai^p V -P-e 



-1/2 



(8) 



Results are collected in Table [V] for a scattering angle of 
5 degrees and for a running time of T = 30 days. At 
E = 1.05 GeV, AR n /R n = 0.66 %. Alternatively, since 
the error scales with T -1 / 2 , a 1% statistical error on R n 
can be obtained in 30 x (0.66) 2 = 13 days. We emphasize 
that this does not include many other sources of error in 
addition to statistics. For example there will be an addi- 
tional normalization error associated with an uncertainty 
AP in the beam polarization P. This will increase the 
total error ARl° l /R n by a factor b from AR n /R n of Eq. 



AR^ = b AR n 



with b, 



b = 



Rn 



■AR„e Rjl ■ 
R n 



(9) 



(10) 



At present P can be determined with an error AP of 
order 1% to 2%. 

This estimate of 13 days for a 1% determination of 
R n may be slightly optimistic compared to the actual 
PREX experiment because the experimental acceptance 
includes a range of scattering angles from 4.5 to 7 de- 
grees and the angle averaged A pv may be slightly less 
sensitive to R n than en n at 5 degrees. Furthermore the 
actual beam current could be less than 100 fiA. Never- 
theless, our estimate of 13 days for a 1% measurement of 
R n in 208 Pb provides a benchmark that is based on our 
assumptions. In Section |IV| we will compare this to the 
running time for some other neutron density measure- 
ments. 

Note that the results in Table [V] for 208 Pb at larger 
energies than 1.05 GeV assume a given shape for the 
neutron density. At higher energies one is also increas- 
ingly sensitive to the surface thickness, see below, and 
other features of the neutron density in addition to the 
radius. Therefore one needs to interpret the sensitivity 
to R n given in Table |V| with care. 



TABLE V: Statistical error estimates for measuring R n in 30 
days. Results are first presented for 208 Pb, 48 Ca, and 40 Ca 
at a laboratory scattering angle of five degrees and then 48 Ca 
results are also presented for a scattering angle of four degrees, 
see text. The neutron and proton densities are calculated in 
the Skyrme HF theory with the SLY4 interaction. 



Nucleus 


E 
GeV 


ppm mb/str 


Rate(5°) e Rrl AR n /R n 
MHz/arm % 


208p b 


1.05 


0.7188 1339 


1736 -2.762 0.6637 


48 Ca 


1.80 


2.358 8.630 


164.3 -4.266 0.4258 


40 Ca 


1.90 


2.301 5.832 


111.0 -3.920 0.5777 




E 




Rate(4°) e fl „(4°) AR n /R n 


48 Ca 


2.20 


2.290 16.56 


315.2 -3.961 0.3409 



IV. ADDITIONAL NEUTRON DENSITY 
MEASUREMENTS 

In this section we consider other parity violating neu- 
tron density measurements. First, in Subsection |IV A| we 
explore a measurement of the surface thickness a n of the 
neutron density in 208 Pb. We assume this measurement 
is after the neutron radius R n has been measured. In 
Subsection |IVB| we look a t R n a nd a n measurements for 

we consider the neutron 
Sn, and 124 Sn. 



Ca, while in Subsection IV C 
radius in Tin isotopes 112 Sn, 12 " 



a Wood Saxon form, 

Pn{r) = Po/[l + cxp(r - Ro)/a n ]. (11) 

An approximate fit of Eq. [TT]to the SLY4 neutron density 
yields a, 

density and the sensitivity, 



0.55 fm. We calculate A pv for this neutron 



din A. 



pv 



j i (12) 
am a n 

to small changes in a n . We assume that the neutron 
r.m.s. radius R n has been fixed by an earlier measure- 
ment. Therefore we calculate the derivative in Eq. [F2|a t 
fixed R n and not at fixed Rq (the parameter in Eq. 111. 
Indeed as a n is changed Rq is also changed so that R n re- 
mains fixed. Figure [5] shows e an for beam energies of 1.05 
and 1.8 GeV. The statistical sensitivity of a measurement 
of a„ is 



= UVi 



total 



A 2 P 2 f 2 

pv c a T 



-1/2 



(13) 



Results for Aa n are collected in Table |VI| In 30 days of 
running at 1.8 GeV and 5 degrees one can obtain a sta- 
tistical error of 7.9 %, using the experimental parameters 
in Table |III| We emphasize that this is a sensitivity to 
a n of 7.9 % given our assumptions that R n is precisely 
known and that the neutron density has a Wood Saxon 
form. 
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FIG. 5: (Color on line) Sensitivity of the parity violating 
asymmetry A pv for 208 Pb to changes in the neutron surface 
thickness a„, see Eq. |11| where e a „ = e a „ = d In A pv jd In a n 
versus scattering angle 6 at beam energies of 1.05 GeV (solid 
line) and 1.8 GeV (dashed line). The SLY4 neutron density 
was used. 



Surface thickness a n in Pb 



Neutron Densities of 48 Ca and 40 Ca 



To explore a surface thickness measurement of the neu- 
tron density in 208 Pb we model the neutron density with 



The neutron density of 48 Ca is very interesting. First 
48 Ca is one of the lightest neutron rich closed shell nuclei, 
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TABLE VI: Statistical error estimates for measuring a n in 
30 days for nuclei 208 Pb, assuming a n w 0.55 fm, and 4S Ca, 
assuming a n w 0.50 fm. The proton density is calculated in 
the Skyrme HF theory with SLY4 and the neutron density is 
of Wood-Saxon form, p(r) = po/[l + ea;p( r ~ fl ° )]. 



Nucleus 


E 
GeV 


A pv (5°) ^(5°) 
ppm mb/str 


Rate(5°) e a „(5°) Aa n /o„ 
MHz/arm % 


208p b 


1.80 


2.085 40.59 


52.61 -0.4582 7.924 


48 Ca 


2.15 


2.744 1.103 


21.00 1.474 2.962 



and is much lighter than 208 Pb. Therefore i?„( 48 Ca) may 
provide independent information that is not contained in 
i? n ( 208 Pb). In addition, more microscopic calculations, 
using for example coupled cluster or no core shell model 
approaches, may soon be feasible for 48 Ca. These are 
presently not feasible for 208 Pb. This may allow one to re- 
late neutron density measurements in 48 Ca more closely 
to microscopic nucleon-nucleon and three nucleon forces. 
In particular, one would like to probe poorly known three 
neutron forces with neutron density measurements. Fur- 
thermore, 48 Ca is a relatively light double beta decay 
nucleus. This further motivates nuclear structure studies 
that could also probe its neutron density. 

Finally, in addition to 48 Ca, the isotope 40 Ca is a sta- 
ble N = Z nucleus where it is expected that R n should 
be close to, but slightly smaller than, R p . Therefore one 
can compare a variety of measurements on the two Ca 
isotopes. Indeed, a parity violating measurement of R n 
for 40 Ca provides a fundamental check on the whole pro- 
cedure to use parity violating electron scattering to mea- 
sure neutron densities. Many systematic errors will can- 
cel in combining two parity violating measurements to 
determine i?„( 48 Ca) — i?„( 40 Ca). For example, radiative 
corrections are expected to be very similar for the two 
isotopes. 

The sensitivity of A pv to R n for 48 Ca is shown in Fig. 
[6j At E = 1.8 GeV and 8 = 5 degrees one can measure 
R n with a statistical sensitivity of 0.43 % in 30 days, 
see Table |Vj Alternatively, one can measure R n with a 
statistical error of 1% in only 5.5 days. This is less than 
half of the beam time required for a 1% measurement 
in 208 Pb, as we discuss below. Table |V| also shows that 
R n in 40 Ca can be measured at E = 1.90 GeV and 8 = 5 
degrees to 0.6% in 30 days, or a 1% error is possible in 11 
days. Calcium 40 is slightly harder to measure than 48 Ca 
because A pv is smaller, given the smaller weak charge, 
and because the sensitivity en n is slightly smaller. 

It is interesting to compare 48 Ca to 208 Pb. First, Table 
|V| shows that the optimal energy, and momentum trans- 
fer q, for a measurement in 48 Ca is higher than in Pb. 
This immediately follows because R n in Ca is smaller 
than lead so the product qR n stays approximately con- 
stant. This higher q insures that A pv is larger for the 
Ca measurement, since A pv scales approximately with 
q 2 . Finally the Ca figure of merit is larger than that for 
Pb because of the higher A pv and because en n is slightly 



larger. 

All of the above results assume a laboratory scattering 
angle 8 = 5 degrees. The existing PREX septum bends 
1.05 GeV electrons, scattered at 5 degrees, into the High 
Resolution Spectrometers. Note that this septum may 
not have a strong enough magnetic field to work at higher 
energies. With the planned energy upgrade at Jefferson 
Laboratory, 2.2 GeV may be a good energy for a neutron 
density measurement. This would allow a single pass 
beam to be used for the measurement while, at the same 
time, higher energy multi-pass beams are sent to other 
experiments. A 2.2 GeV neutron radius measurement on 
208 Pb may need detectors at very small angles just over 
two degrees. However a 2.2 GeV R n measurement on 
48 Ca only needs detectors near four degrees. Therefore 
Table [v] also lists results for 48 Ca at four degrees. Be- 
cause of the high beam energy the figure of merit is very 
good. In 30 days one could, in principle, reach a sensitiv- 
ity of 0.34 %. Alternatively, one only needs 3.5 days to 
get 1 % statistics for R n . This is a factor of four shorter 
time than for PREX. If a new septum could be designed 
to work at four degrees and 2.2 GeV, this would allow 
a large improvement in the figure of merit over PREX. 
Furthermore, many systematic errors depend on the ab- 
solute size of A pv . These may be easier to deal with 
in a 48 Ca measurement, than for PREX, because A pv is 
larger. 




© (deg) 

FIG. 6: (Color on line) Sensitivity of the parity violating 
asymmetry A pv for 48 Ca to changes in the neutron radius 
ER n = €R n = din A pv /din R n versus scattering angle 6 at 
beam energies of 1.8 GeV (solid line) and 2.15 GeV (dashed 
line). The SLY4 neutron density was used. 

The sensitivity of A^ v to the surface thickness a n for 
48 Ca is shown in Fig. g At E = 2.15 GeV and 6 = 5 de- 
grees, one is sensitive to a n with a statistical error of 3.0% 
after 30 days, see Table |VI| This is much smaller than 
the 7.9% error for a n in 21)8 Pb. There are two reasons 
for this dramatically increased sensitivity. First, A pv is 



large, 2.7 ppm for 48 Ca at 2.15 GeV. Second, one is much 
more sensitive to a n in 48 Ca than in 208 Pb because 48 Ca 
is mostly surface while the surface is only a small part of 
208 Pb. This leads to a much larger e an . 
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FIG. 7: (Color on line) Sensitivity of the parity violating 
asymmetry A pv for 4 8 Ca to changes in the neutron surface 
thickness a n , see Eq. 11 e a „ = e<j„ = dm Aj,„/dlna n versus 
scattering angle 8 at beam energies of 1.8 GeV (solid line) 
and 2.15 GeV (dashed line). 



Finally Table VIII presents statistical error estimates for 
measuring the surface thickness a n in 120 Sn. This could 
be measured at E = 2.15 GeV and 9 = 5 degrees with a 
statistical sensitivity of 5.1% after 30 days. Again this is 
a smaller error than for 208 Pb. Finally, we expect sim- 
ilar results for an a n measurement in other even N Sn 
isotopes. 



TABLE VII: Statistical error estimates for measuring R n for 
112 Sn, 120 Sn, and 124 Sn in 30 days. The neutron and proton 
densities are calculated in the Skyrme HF theory with the 
SLY4 interaction. 



Nucleus 


E 
GeV 


A pv (5°) &(5°) 
ppm mb/str 


Rate(5°) e Rn (5°) AR n /R n 
MHz/arm % 


112 Sn 


1.30 


1.099 187.2 


654.7 -3.157 0.6183 


120 Sn 


1.25 


1.124 230.9 


807.8 -3.070 0.5599 


124 Sn 


1.25 


1.160 223.4 


781.4 -3.172 0.5337 



TABLE VIII: Statistical error estimates for measuring a n 
in 30 days for 120 Sn, assuming a n m 0.55 fm. The proton 
density is calculated in the Skyrme HF theory with SLY4 
and the neutron density is of Wood-Saxon form, p(r) = 
po/[l + exp{ 1 



-Ra * 



Nucleus 


E 


A pv (5°) 


It(5°) 


Rate(5°) e a „(5°) Aa n /a n 




GeV 


ppm 


mb/str 


MHz/arm % 


120 Sn 


2.15 


3.044 


4.470 


15.63 -0.8889 5.131 



C. Neutron Densities of 112 Sn, 120 Sn, and 124 Sn 

The neutron density of Tin isotopes are interesting for 
a number of reasons. First a variety of Tin isotopes are 
available for experiments with either stable or radioac- 
tive beams. Indeed heavy ion collisions with Tin isotopes 
have been used to probe the density dependence of the 
symmetry energy |31| . The neutron radius of a heavy nu- 
cleus has been shown to be closely related to the density 
dependence of the symmetry energy. A rapid density de- 
pendence implies a large pressure for neutron rich matter 
and this give a large neutron radius [23], Second, pairing 
corrections may be larger in Tin isotopes than in a closed 
shell nucleus such as 208 Pb, and these corrections could 
impact R n . 



Tabic VII presents statistical error estimates for R n 
112 Sn, 



measurements in 



J Sn, and Sn. In general 
there is a smooth dependence on neutron number TV with 
larger N isotopes being somewhat easier to measure be- 
cause they have larger weak charges and A pv . We expect 
similar results for other even N Sn isotopes. We see 
that R n can be measured in 120 Sn at E = 1.25 GeV and 
8 = 5 degrees with a statistical sensitivity of 0.56 % after 
30 days. This is a slightly smaller error than for 208 Pb. 
Although 120 Sn has a smaller cross section than 208 Pb, 
the asymmetry A pv is larger for 120 Sn than for 208 Pb and 
this more than compensates for the smaller cross section. 



V. DISCUSSION AND CONCLUSIONS 

Parity violating electron scattering allows model inde- 
pendent measurements of neutron densities that are free 
from most strong interaction uncertainties. In this pa- 
per we present statistical error estimates for a variety of 
experiments. The neutron radius R n can be measured 
in several nuclei, as long as the nuclear excited states 
are not too low in energy. In general, we find that the 
smaller the nucleus, the easier the measurement. This is 
because smaller nuclei can be measured at higher momen- 
tum transfers where the parity violating asymmetry A pv 
is larger. Also in general, the more neutron rich the iso- 
tope, the easier the measurement, because neutron rich 
isotopes have larger weak charges and larger A pv . 

Since measurements of R n are feasible in many nu- 
clei, one can use this freedom to choose an element that 
makes a very robust target. Alternatively, PREX uses 
208 Pb because of its very simple nuclear structure. The 
doubly magic 208 Pb is an excellent closed shell nucleus 
where a variety of corrections, such as pairing correla- 
tions, may be small. This allows a clean interpretation 
of a R n measurement in terms of bulk properties of neu- 
tron rich matter such as the equation of state [12"] . 

However, in this paper, we find that measuring R n in 
48 Ca is very interesting because 48 Ca is much smaller 
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than Pb. As a result R n can be measured faster with 
a higher figure of merit. In addition, i?„( 48 Ca) may 
be more easily related to two nucleon and three nu- 
cleon interactions, including very interesting three neu- 
tron forces, than i?„( 208 Pb). This is because 48 Ca has 
fewer nucleons than 208 Pb and this greatly simplifies mi- 
croscopic coupled cluster or no core shell model calcu- 
lations. In considering a parity violating electron scat- 
tering experiment on 48 Ca, one should also consider the 
information that may already be available on the distri- 
bution of f 7/2 neutrons from inelastic magnetic electron 
scattering. However the interpretation of this inelastic 
data may be somewhat model dependent because of nu- 
clear structure effects such as core polarization. 

After measuring R n , one can constrain the surface 
thickness of the neutron density a n with a second mea- 
surement at somewhat higher momentum transfers. We 
present error estimates for measuring a n in 48 Ca, 120 Sn, 
and 208 Pb. Again, we find that a n is easier to measure 
in smaller nuclei. One should study the sensitivity of 
the surface thickness to different features of the effective 
interaction to determine the possible nuclear structure 
information that would be available from a measurement 
of the surface thickness. 

Finally in future work, we will present statistical er- 
ror estimates for using measurements at several momen- 
tum transfers to determine the complete neutron density 
p n (r) in a model independent fashion. This appears fea- 
sible, but difficult, for 48 Ca. However determining p n (r) 
may be extremely difficult for 208 Pb. This measured neu- 
tron density for 48 Ca, combined with the previously mea- 
sured charge density, will provide a very detailed picture 



of an atomic nucleus. 

We find that parity violation experiments are feasible 
for a variety of nuclei and that experiments on lighter 
nuclei have larger figures of merit, in general. This fur- 
ther motivates studies of correlations between the neu- 
tron radii of a variety of nuclei, see for example [32j . 
These studies should consider many nonrelativistic and 
relativistic effective interactions and include deformation, 
pairing, and other nuclear structure effects. If these stud- 
ies show that neutron radii in heavy and medium light 
nuclei are not strongly correlated than an additional par- 
ity experiment in a lighter nucleus such as 48 Ca should 
provide new nuclear structure information in addition to 
that provided by PREX for 208 Pb. 

Alternatively, if neutron radii in heavy and light nu- 
clei are strongly correlated, than our results may be even 
more important. In this case one may have considerable 
freedom to choose the nucleus to minimize the exper- 
imental difficulties of a measurement. For example, by 
following up the PREX experiment with a 48 Ca measure- 
ment one may be able to significantly improve the accu- 
racy of the determination of the neutron radius, beyond 
what is actually achieved in PREX. This is because the 
larger figure of merit allows more statistics to be accumu- 
lated in a shorter time, while the larger parity violating 
asymmetry may reduce some systematic errors. A more 
accurate determination of the neutron radius could be 
very useful to better constrain effective interactions for 
nuclear structure. 

This work was supported in part by DOE grant DE- 
FG02-87ER40365. 
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